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We present partial branching fractions for inclusive charmless semileptonic B decays B — * X u lv, 
and the determination of the CKM matrix element \V u b\- The analysis is based on a sample of 
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383 million T(4S) decays into BB pairs collected with the BABAR detector at the PEP-II e + e~ 
storage rings. We select events using either the invariant mass Mx of the hadronic system, the 
invariant mass squared, q 2 , of the lepton and neutrino pair, the kinematic variable P+ or one of 
their combinations. We then determine partial branching fractions in limited regions of phase space: 
AB= (1.18±0.09 st at. ±0.07 syst . ±0.01theo.) x 1(T 3 {M x < 1.55 GeV/c 2 ), AB = (0.95 ± O.lfW ± 
0.08 syst .±0.01th e o.)xl0- 3 (P+ < 0.66 GeV/c),andAS = (0.81±0.08 st at.±0.07 syst .±0.02 thco .)xl0- 3 
(M x < 1.7 GeV/c 2 , q 2 > 8 GeV 2 /c 4 ). Corresponding values of \V ub \ are extracted using several 
theoretical calculations. 

PACS numbers: 13.20.He, 12.15.Hh, 14.40.Nd 



In the Standard Model the element V u b of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [1] 
plays a critical role in tests of the prediction of CP vi- 
olation. Since the rate for charmless scmilcptonic de- 
cays, B — > X u lv [2], is proportional to |Vu(,| 2 , and the 
hadronic and leptonic currents are factorizable, the best 
method to extract this quantity is to measure branching 
fractions for such decays [3]. Experimentally, the princi- 
pal challenge is to separate the rare B — > X u lv decays 
from the approximately 50 times larger B — > X c lv back- 
ground. Given that the u quark is much lighter than the 
c quark, regions of phase space can be defined where the 
background is suppressed. To relate the decay rate of 
the B meson to \V u b\ 1 parton level calculations have to 
be corrected for perturbative and non-perturbative QCD 
effects. A variety of QCD calculations are available to 
determine these corrections [4-6]. 

In this letter, we present a measurement of partial 
branching fractions for inclusive charmless scmilcptonic 
decays, B — > X u lv [7]. T{AS) — > BB events are tagged 
by the full reconstruction of a hadronic decay of one of 
the B mesons (J5 reco ). The semileptonic decay of the sec- 
ond B meson (B rC coii) is identified by the presence of an 
electron or a muon. This technique results in a low event 
selection efficiency but allows the determination of the 
momentum, charge, and flavor of the B mesons. 

We use three kinematic variables to separate B —* 
X u lv decays from the dominant B — > X c lT> background: 
Mx, the invariant mass of the hadronic system X u , c ; 
q 2 , the invariant mass squared of the lepton-neutrino 
system; and P + = Ex — \Px\ [4, 5], where Ex and 
Px are the energy and momentum of the hadronic sys- 
tem X UtC calculated in the B rest frame. We measure 
the fraction of partial rates of charmless semileptonic 
decays Ai? u/sl = AB(B -> X U £9)/B(B -> Xtv) in re- 
stricted phase space regions, corrected for resolution ef- 
fects. The resulting partial branching fractions are used 
to calculate \V u b\ following theoretical prescriptions. 

The analysis uses a sample of 383 million T(AS) de- 
cays into BB pairs, corresponding to an integrated lu- 
minosity of 347.4 fb _1 , collected with the BABAR de- 
tector [8]. Charmless semileptonic B — > X u lv decays 
are simulated as a combination of three-body decays 
(X u = 7T, 77,77' , p,ui, . . .) [9] and decays to non-resonant 
hadronic final states X u [10]. The motion of the b quark 



inside the B meson is modeled with the shape function 
parametrization given in Ref. [10]. The simulation of the 
B — ► XJLv background uses an HQET parametrization 
of form factors for B — > D*tv [11, 12], and models for 
B -» Dirtv, D*ir£v [13], and for B -> D£U,D**£p [9]. 
The simulation of the hadronization is performed by 
Jetset7.4 [14]. We use GEANT4 [15] to simulate the 
detector response. 

To reconstruct a large sample of hadronically decaying 
B mesons, B rcco — ► £)(*)y ± are selected. Here, the sys- 
tem Y ± consists of hadrons with a total charge of ±1, 
composed of ni7r ± n 2 K ± n 3 K® n 4 7r°, where m + n 2 < 5, 
n 3 < 2, and 714 < 2. The kinematic consistency of 
B reC o candidates is checked with two variables, toes = 
y/s/A. — p% and AE = E B — y/s/2. Here yfs is the total 
energy in the T(AS) center of mass frame, and ps and 
Eb denote the momentum and energy of the -Brcco candi- 
date in the same frame. We require AE — within three 
standard deviations as measured for each decay mode. 
For each of the -B reC o decay modes, the purity V is es- 
timated using Monte Carlo (MC) simulation. V is de- 
fined as the ratio of signal over background events with 
n^ES > 5.27GeV/c 2 . Only modes for which V exceeds 
20% are used. On average, we reconstruct at least one B 
candidate in 0.3% (0.5%) of the B°B° (B + B~) events. 
For events with more than one reconstructed B decay, 
the decay mode with the highest purity is selected. 

We determine the number of -B r cco candidates from an 
unbinned maximum likelihood fit to the toes distribu- 
tion. The data are fit to the sum of three contributions: 
signal Brcco decays, combinatorial background from BB 
events, and continuum (e+e~ — > qq, q = u, d, s, c) events. 
A Threshold function [16] is used to describe the combi- 
natorial and continuum backgrounds. To obtain a good 
description of the signal mEs distribution, we adopt the 
modified Gaussian function used in Ref. [17], to account 
for energy losses of photons in the detector. Fits to the 
toes distribution are shown in Fig. 1. Semileptonic de- 
cays B — ► Xiv of the Brecoii candidate are identified by 
an electron or muon with momentum, p\, defined in the 
B rest frame, greater than 1 GeV/c. For charged _B r cco 
candidates, we require the charge of the lepton to be con- 
sistent with a prompt semileptonic B decay. For neutral 
Brcco candidates, both charge-flavor combinations are re- 
tained and the known average B°-B° mixing rate [18] is 
used to extract the prompt lepton yield. 
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FIG. 1: The to.es distribution for data (full circles) is shown 
together with the results of the fit (solid line) for selected 
semileptonic decays from B + B~ events (a) and B°B° events 
(b). The dashed line shows the contribution from combinato- 
rial and continuum background. 



The hadronic system X in the decay B — ► Xlv is 
reconstructed from charged tracks and energy deposi- 
tions in the calorimeter that are not associated with 
the -B roco candidate or the identified lcpton. We re- 
construct by performing a mass-constrained fit to 
7r + 7r~ pairs with an invariant mass in the range 0.473- 
0.523 GeV/c 2 . The neutrino four-momentum p v is esti- 
mated from the missing momentum four-vector p m i ss = 
Pr(4S) — PB teco ~Px~Pe, where all momenta are measured 
in the laboratory frame and Pr(4S) refers to the T(AS) 
meson. 

To select B — > X u lv candidates we require exactly one 
charged lepton with p\ > lGeV/c, charge conservation 
(Qx + Qi + QB roC o = 0); an d a missing mass consistent 
with zero {m 2 niss < 0.5GeV 2 /c 4 ). These criteria sup- 
press the dominant B — > X c lv decays, many of which 
contain additional leptons or an undetected K® meson. 
We suppress the B — > D*lv background by reconstruct- 
ing the low momentum 7r + from the D* + — > D°7r + decay. 
Since the momentum of the 7r + is almost collincar with 
the D* + momentum p^*+ , we can approximate the D* + 
energy as E D »+ ~ m D »+ x £' 7r /145 McV/c 2 . The neutrino 
mass TO 2 cto = (pb—Pd-+ ~Pe) 2 ls peaked at zero for back- 
ground events. The requirement TO 2 eto < — 3GeV 2 /c 4 
reduces the B — > background by about 36% while 

keeping more than 90% of signal events. We reject events 
with charged kaons or in the -B r0C oii to reduce the 
background from B — ► decays. 

To extract the distribution in the variables Mx, -P+, 
and the combination of Mx and q 2 , we perform fits 
to the -Brcco toes distributions for subsamples of events 
in individual bins for each of the variables, and subse- 
quently separate the signal from the combinatorial and 
continuum backgrounds for the three distributions. The 
resulting distributions are presented in Fig. 2. To re- 
duce the systematic uncertainties in the derivation of the 
branching fractions we determine the ratios of the par- 
tial branching fractions to the total semileptonic branch- 



ing fraction. This is done for restricted regions of phase 
space, M x < 1.55 GeV/c 2 , P + < 0.66 GeV/c, and 
(M x < 1.7 GeV/c 2 , q 2 > 8.0 GeV 2 /c 4 ). Specifically 
we define this ratio as 



AB(X u £De) (N u - N™* - BG u )/(e« el e£ n ) 



B(Xlv t ) 



(N sl - BG A ) 



si si 



(1) 

where N u refers to the number of observed events, BG U 
to the estimated number of background events, and iV° ut 
to the signal events that migrate from outside the kine- 
matic region into the signal region. They are deter- 
mined by a x 2 fit to the measured spectra with signal 
and background shapes determined from MC simulation. 
N sl = 181074 ± 706 and BG sl = 12185 ± 78 are the 
number of semileptonic events, extracted with a mgg 
fit, and the corresponding background, determined from 
simulation. The efficiency e" el denotes the fraction of se- 
lected -B r eco-tagged signal events with a high-energy lep- 
ton. The model-dependent efficiency ej^ in accounts for the 
loss of selected events generated in the kinematic region 
that migrate outside this region. The efficiency of the 
tag and lcpton selection, e t and q, differ slightly for the 
signal and the semileptonic samples, due to differences in 
the lepton momentum distribution and the multiplicity 
of the recoiling B meson. To convert the ratio in Eq. 1 to 
partial branching fractions, we use the total semileptonic 
branching fraction, B(B -> XlV e ) = (10.75±0.15)% [18]. 
The resulting partial branching fractions for the three se- 
lected kinematic regions, along with parameters in Eq. 1, 
are listed in Table I. The statistical correlations between 
the M x and {M x ,q 2 ), P+ analyses are 65%, 67%, 38% 
respectively. 

We consider several sources of systematic uncer- 
tainties. Detector- related uncertainties take into ac- 
count particle (e, fj,, K) identification (efficiency, mis- 
identification), charged particle tracking efficiency, pho- 
ton reconstruction efficiency and K® interactions. We 
estimate the uncertainty due to signal and background 
modeling. The uncertainty on the signal modeling are 
due to the modeling of exclusive charmless semileptonic 
decays and gluon splitting into ss-quark pairs. We also 
calculate the uncertainties due to the non-perturbative 
parameters and the functional form of the shape func- 
tion. The background simulation depends on the B and 
D branching fractions and B — > D*lv form factors; the 
corresponding systematic uncertainties are calculated by 
varying all these quantities within their experimental er- 
rors. We estimate the error due to toes fits, coming from 
the uncertainty in the parameterization ansatz. Finally, 
we estimate the error due to MC statistics. The fractional 
contribution of each uncertainty is shown in Table II to- 
gether with the total error. 

The results of the partial branching fractions are trans- 
lated into \V u b\ in the context of recent QCD calcula- 
tions [4-6] , including estimates of theoretical uncertain- 
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FIG. 2: Upper row: measured Mx (a), P+ (b) and q 2 with Mx < 1.7 GeV/c 2 (c) spectra (data points). The result of the fit 
to the sum of three MC contributions is shown in the histograms: B — > X u tv decays generated inside (no shading) and outside 
(dark shading) the selected kinematic region, and B — > X c tv and other background (light shading). Lower row: corresponding 
spectra for B X u tv after B — ► X c lv and other background subtraction; they have been rebinned in order to show the shape 
of the kinematic variables. 



TABLE I: Summary of the fitted number of events and efficiencies, AB(B — ► X u £v), and extracted \V u b\ f° r the three kinematic 
cuts. The first uncertainty is statistical, the second systematic. For AB, the third uncertainty is due to the theoretical 
knowledge of the signal efficiency; for the \V n b\ values, it comes from the the theoretical uncertainty on A(. For Ref. [4] we use 
the exponential parametrization of the shape function. 



Method 


N u 


AfOUt 


BG U 


U U 

scl kin 


' 1 ' t 


AB(B -> X u lv) (10" 3 ) 


\V ub \ x (10" 3 ) 


M x 


803 ± 60 


27 ± 2 


923 ± 21 


0.331 ±0.003 


0.76 ±0.02 


1.18 ±0.09 ±0.07 ±0.01 


4.27 ±0.16 ±0.13 ±0.30 [4] 
4.56 ±0.17 ±0.14 ±0.32 [5] 


P+ 


633 ± 63 


48 ±5 


1183 ± 27 


0.344 ± 0.003 


0.81 ±0.02 


0.95 ±0.10 ±0.08 ±0.01 


3.88 ±0.19 ±0.16 ±0.28 [4] 
3.99 ±0.20 ±0.16 ±0.24 [5] 
















4.57 ±0.22 ±0.19 ±0.30 [4] 


M x ,q 2 


562 ± 55 


32 ± 2 


789 ±9 


0.353 ± 0.005 


0.79 ±0.03 


0.81 ±0.08 ±0.07 ±0.02 


4.64 ±0.23 ±0.19 ±0.25 [5] 
4.93 ± 0.24 ± 0.20 ± 0.36 [6] 



ties (see Table I). The hadronic input parameters, the 
6-quark mass mb, and the kinetic energy expectation 
value /it 2 , are extracted from moment measurements in 
B — > X s j and B — > X c iv. Their values in the ki- 
netic scheme [19] are m b = (4.59 ± 0.04) GeV/c 2 and 
fil = (0.40 ± 0.04) GeV 2 /c 2 [20] and are translated into 
values in different schemes, as needed [4-6]. The partial 
branching fraction AB(B — > X u £p) is related directly to 
\V ub \ by the relation \V ub \ = [AB(B -> XJ^/nAC} 1 / 2 , 
where r b is the average B lifetime [18], and AC is the 
prediction for the partial rate for B — » X u tv in the given 
phase-space region [4-6]. 

In summary, we have measured the branching fractions 
for inclusive charmless semileptonic B decays B — > X u iv 
in three overlapping regions of phase space. Relying on 



theoretical predictions, we extract values for the CKM 
matrix element \V ub \ from our measured AB. 

We find that the determinations of \V ub \ agree at 1 a 
level in the BNLP framework for the Mx and combined 
(Mx,q 2 ) analyses. The analysis based on P + differs from 
the two others at a 2.5 a level, as indicated also by other 
experiments [21]. The Mx analysis captures the largest 
portion of phase space and gives the most precise deter- 
mination of \V u b\- Within their stated theoretical uncer- 
tainties, the results based on BLNP and DGE give con- 
sistent results. The result, based on the hadronic mass 
spectrum, supersedes our previously published measure- 
ment [3], reducing the relative uncertainty by 40%. These 
values are in good agreement with other inclusive IV^&I 
determinations and they are somewhat higher, though 
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TABLE II: Contributions to the systematic uncertainty on the measured AB(B — > X u lv), shown in percent (%) for the three 
kinematic cuts, from: detector, shape function (input parameters and functional form), exclusive B(B — ► X u tv), gluon splitting, 
exclusive B(B — > X c iv), B — ► D*£~V form factors, B(D), tties fit, MC statistics. The last column gives the total systematic 
uncertainty. 



Method Detector sha P e 
function 


B(B 


= 7T,p, . . . 


Gluon 
splitting 


B(B -> A C £P) 


B -» 

form factors 


B(£>) 


m B s fit Mo f nt f ( ? 9 arl ° Total 
statistics 


M x 1.92 0.90 




2.08 


1.62 


0.87 


0.21 


0.44 


3.71 3.22 6.07 


P+ 3.88 1.31 




2.22 


1.47 


2.80 


0.39 


0.73 


3.98 4.62 8.38 


M x , q 2 3.83 2.43 




2.71 


1.02 


1.17 


0.55 


0.79 


5.17 4.29 8.81 



compatible, than the results based on exclusive charm- 
less semileptonic decays [18]. 
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